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ABSTRACT 

“The Duck” is a complicated non-thermal radio system, consisting of the energetic radio pulsar B1757- 
24, its surrounding pulsar wind nebula G5.27-0.90 and the adjacent supernova remnant (SNR) G5.4-1.2. 
PSR B1757-24 was originally claimed to be a young (w 15000 yr) and extreme velocity (> 1500 km s“') 
pulsar which had penetrated and emerged from the shell of the associated SNR G5.4-1.2, but recent upper 
limits on the pulsar’s motion have raised serious difficulties with this interpretation. We here present 8.5 GHz 
interferometric observations of the nebula G5.27-0.90 over a 12-year baseline, doubling the time-span of pre¬ 
vious measurements. These data correspondingly allow us to halve the previous upper limit on the nebula’s 
westward motion to 14 milliarcseconds yr“' (5 -(t), allowing a substantive reevaluation of this puzzling ob¬ 
ject. We rule out the possibility that the pulsar and SNR were formed from a common supernova explosion 
~ 15000 yrs ago as implied by the pulsar’s characteristic age, but conclude that an old (> 70000 yr) pulsar 
/ SNR association, or a situation in which the pulsar and SNR are physically unrelated, are both still viable 
explanations. 

Subject headings: ISM: individual (G5.4-1.2) — pulsars: individual (B1757-24) — radio continuum: ISM — 
stars: neutron — supernova remnants 


1. INTRODUCTION 


The formation and subsequent evolution of pulsars are not 
yet fully understood. A powerful constraint on these pro¬ 
cesses is provided by an independent age estimate. In cases 
where both a pulsar and its associated supernova remnant 
(SNR) can be identified, an age estimate which is independent 
of both distance and inclination effects is simply tp = 0//i, 
where 0 is the angular separation between the pulsar and 
its inferred birth site, and ix is the pulsar ’s proper motion 
llMigliazzo et al.l 120021 iKramer et al.ll200^ . This can then 
be compar ed to the age of the system as expected from spin- 
down te.g.. lUorimer & KrameTi200.5h : 


P 

(n-l)P 



( 1 ) 


where P is the current spin period, P is the time-derivative 
of P, P() is the period at birth, and n is the “braking index” 
(see further discussion in 33.2> . Comparison of these two in¬ 
dependent age estimates provides information on Pq and n, 
i.e., on the processes which impart and then dissipate the con¬ 
siderable angular momentum of neutron stars (e.g., Gaensler 
& Erail 2000, hereafter GEOO; Kaspi et al. 2001b). If one 
assumes that Pq ^P and n = 3, Equation m reduces to the 
expression for the “characteristic age” of a pulsar, Tc =Pj2P. 
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PSR B1757-24 is an isolated 125-ms pulsar, surrounded 
by the cometary radio and X-ray pulsar wind nebula (PWN) 
G5.27-0.90, w hich in turn is located just outside the 
SNR G5.4-] .2 jCaswell et al Jll 987t iManchester et allfT^ 

iFrail & Kulkarnil199lHFrail et a1.l1994hHKasr)i et a1.l2nnil) . 

As shown in Figure ^ the combined system has a distinc¬ 
tive morphology which has led to it being termed “the Duck”. 
Because of the proximity of the pulsar to the SNR, and be¬ 
cause the PWN morphology suggests that the pulsar is mov¬ 
ing away from the SNR interior, it has been widely assumed 
that the pulsar and the SNR are physically associated. The 
SNR and PWN also have very similar HI absorption spe ctra, 
both suggesting a distance of ^ 5 kpc JFrail et al.ll994bl) . and 
consistent with the distance of 5.1 kpc implied by the pulsar’s 
di spersion measure and t he Galactic electron density model 
of ICordes & Laziolj2002ft . We consequently adopt a common 
distance of 5 kpc for pulsar, PWN and SNR in further discus¬ 
sion. 

PSR B1757-24 has a characteristic age = 15 kyr 
JManchester et alJl2005l) . If we assume tp = Tc and that the 
pulsar has travele d 0 = 16fl-20.'6 from the SNR’s cen¬ 
ter in its lifetime JFrail et alJll994bl) .^ this implies a west¬ 
ward proper motion for the pulsar of magnitude ^ = 0 /tc = 
63 - 80 mas yr“^. However, radio interferometric observa¬ 
tions of the western tip of the PWN taken at the Very Large 
Array (VLA) over 6.7 years yielded a surprising 5-cr upper 
limit fj, < 25 mas yr“' (GEOO), implying a pulsar age tp > 
39-50 kyr 3> t^, if the pulsar was born near the SNR’s geo¬ 
metric center. GEOO used the stand-off distance of the PWN, 
the radius of the SNR and the separation of the pulsar from the 
SNR’s center to derive a solution for the system’s evolution 
which predicted ^ 90- 170 kyr and p, 10 mas yr“*. Sub¬ 
sequently, Thorsett et al. (2002, hereafter TBG02) observed 

® The range of values quoted for 0 reflects the uncertainty in locating 
the pulsar’s presumed birthplace. This results from the fact that the pulsar’s 
inferred trajectory does not p ass through the SNR’s geometric center (Fig.m 
IFrail et alJ1994hlt . See ^3.2i for further discussion. 
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the pulsar itself with the VLA over a 3.9-year baseline, and 
placed an independent upper limit on westward proper mo- 
tion'° of /r < 37 mas yr“’. TBG02 argue that this is most 
easily explained if PSR B1757-24 and SNR G5.4-1.2 are un¬ 
related, and if the pulsar is instead moving away from the 
center of the PWN G5.27-0.90. Assuming tp « t^, one can 
then predic t a proper motion u ~ 5 mas yr“'. As a further 
alternative, iGvaramadzd J2004I) has argued that this system 
results from a massive high-velocity progenitor star which 
went supernova inside its wind-blown bubble. The pulsar be¬ 
gan its life substantially offset from the cavity’s center, but 
the resulting SNR e xpands to take on the shape of the cavity. 
IGvaramadzd d2004l) subsequently develops a model in which 
tp sa 54 kyr, 0 « 6' and /r « 7 mas yr“'. 

As an attempt to resolve this puzzling situation, we have 
conducted a new observation of the western cometary tip of 
PWN G5.27-0.90 (see inset to Pig.Q, which doubles the time 
baseline considered by GEOO to 12 years. In ^we present our 
new observations and corresponding measurement of proper 
motion, while in ^we we interpret these new results in the 
context of various possibilities proposed for the origin and 
evolution of this system. 

2. OBSERVATIONS AND ANALYSIS 

We have utilized three separate observations of G5.27-0.90, 
all carried out near 8.5 GHz with the VLA in its hybrid BnA 
configuration. The first two epochs were on 1993 Eeb 02 and 
1999 Oct 23, as discussed by GEOO; the third epoch was on 
2005 Jan 22, with an on-source integration time of 2.9 hr, 
and parameters otherwise the same as for the 1999 epoch 
reported by GEOO. In particular, observations of G5.27- 
0.90 at all three epochs were phase-referenced to the source 
PMN J1751-2524, located 2.°1 from G5.27-0.90 (and for 
consistency using the same position for PMN J1751-2524 
at each epoch, even though the best estimate of the position 
of this source has been updated by a small amount in recent 
years). 

Analysis was carr ied out in the MIRIAD package 
JSault & Killeenll2004l) . Data from the three epochs were re¬ 
duced in almost identical fashion, making allowances for the 
slightly different correlator mode used in the 1993 data. The 
data were edited, calibrated, and imaged using square pixels 
of size 50 mas x 50 mas. The fields were then deconvolved 
using CLEAN and smoothed to a common resolution of 0."85. 
The resulting images, shown in Eigure|2l suggest that there 
have been slight structural changes between epochs. 

The vertical dashed lines in Eigure|2] demonstrate that mo¬ 
tion at the 5-(j limit of GEOO could now have been easily 
discerned, and that any change in the position of the leading 
edge of the PWN is well below this level. We have quantified 
this result by measuring the shift between epochs using the 
MIRIAD task IMDIEE IMDIEE finds the shift that minimizes 
the RMS of intensity fluctuations in the resulting difference 
map, employing cubic convolution i nterpolation t o calculate 
shifts of a non-integer pixel number dPowe 11I1964II . 

Using the approach described by GEOO, the shift that we 
determine between the 1993 and 2005 epochs is 71 ±24 mas 
in a westward direction. However, this does not take into ac¬ 
count any systematic errors present in this shift determination, 
nor does it use the combination of all three epochs. We have 

TBG02 state that their limit is /r < 16 mas yr“* at 95% confidence (i.e., 
2 (t). Here we use 5-cr limits throughout, and have adjusted their limit accord¬ 
ingly. 


thus determined the motion of the PWN by measuring the 
shift between all three possible pairs of epochs using IMD- 
lEE. To characterize the statistical error in the reported shifts, 
for each pair a series of phase shifts were applied to the u-v 
visibility of one epoch, resulting in images which had been 
shifted east or west by 0, 0.4, 0.8, 1.2, ..., 4.0 pixels. This 
allowed us to probe the systematic errors introduced by IMD- 
lEE for non-integer pixel shifts of a faint extended source in 
the presence of noise. To minimize the effects of any internal 
structural change in the PWN between epochs, the brightness 
distribution above a surface brightness of 70 mJy beam“^ was 
pegged at 70 mJy beam"', giving the nebula a largely uni¬ 
form intensity out to its periphery. Eor each pair of epochs 
and each phase shift, IMDIEE was then run using 16 different 
sets of inputs, produced by alternating between four different 
spatial windows, by switching between the choice of epoch 
used as the reference image, and constraining IMDIEE to And 
a shift only in R.A., or in both R.A. and Deck The ensemble 
of outputs from the combinations of these options allowed us 
to characterize, for each input trial shift, a mean fitted shift 
in the R.A. direction and its systematic error. The recovered 
shifts were then plotted against the input shifts and fit with a 
weighted linear least-squares model. The ordinate intercept 
then represented the best-fit measured motion between each 
pair of epochs. The statistical error on this estimate was deter¬ 
mined by calculating the standard error in the mean amongst 
the residuals between input and mean output values for the 21 
trial shifts for each pair. The statistical and mean systematic 
errors were then combined in quadrature to give the final error 
estimate for each pair. Other sources of error, such as phase 
transferral and calibration error, are negligible compared to 
these effects and have not been incorporated in this analysis. 

We then examined the correlation between shifts between 
each of the three epoch pairs and the time separation be¬ 
tween these pairs. If steady westward motion was detected, 
we would have expected these two quantities to correlate, but 
the data show no such trend. We conclude that despite our 
best efforts to minimize the effects of structural changes in 
the nebula between epochs, the underlying proper motion is 
too small in magnitude to be seen in the presence of the sys¬ 
tematic errors introduced by these changes. We have corre¬ 
spondingly determined an upper limit on the source’s motion 
by calculating the of the best linear At to the data, and then 
finding the larger of the two slopes for which increases by 
5-a over this best-fit value. 

Through this approach, we And an upper limit on westward 
proper motion of G5.27-0.90 of /r < 13.9 mas yr“'. This and 
previous estimates are shown in Table Q it can be seen that 
this new result is nearly a factor of two more constraining 
than existing measurements. Eor 0 = 16(1 -20.'6 as assumed 
above, the 5-a lower limit on the system’s age is = 0//i > 
69-88 kyr (see Tabled. 

Note that because the X-ray and radio morphologies of 
G5.27-0.90 suggest proper motion purely in a westward di¬ 
rection, we have not included any shift in Deck in our Anal 
proper motion estimates. In any case, the derived shift be¬ 
tween each pair of epochs is consistent with zero proper mo¬ 
tion in Deck 

Assuming a distance to the system of 5 kpc, we can in¬ 
fer an upper limit on the projected westward velocity for 
PSR B1757-24 of 340 km s“*. This limit is consistent with 
the observed range of mo tions for other pul sars associated 
with SNRs (see Table 6 of lHobbs et akll2()0^ . and also with 
the overall projected velocity distribution of the young pulsar 
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population JArzoumanian et al.l2002HHobbs et al.l2005l) . 

3. DISCUSSION 

Since we have used the motion of PWN G5.27-0.90 as a 
proxy for that of PSR B1757-24, the two objects must move 
together rigidly for our limit to be applicable to the pulsar. As 
discussed by GFOO, the pulsar could be moving faster than our 
upper limit if the stand-off distance between the pulsar and the 
head of the bow shock were steadily decreasing with time, 
as might be produced by the system encountering a sudden 
increase in density. However, this is unlikely to be the case 
since, if anything, the data in Figure |2l suggest that emission 
from the leading edge of the bow shock has slightly faded in 
brightness as a function of time. If this is naively interpreted 
as corresponding to a reduction in ambient density, then the 
stand-off distance could even be increasing with time, further 
strengthening our upper limit on the pulsar’s motion. In any 
case the direct observations of the pulsar reported by TBG02 
rule out a large disparity between the proper motions of the 
pulsar and its PWN. 

The angular proximity and consistent distance estimates for 
PSR B1757-24 and SNR G5.4-1.2 obviously demand that the 
possibility of a physical association be considered. Further 
evidence strengthening the case for a genuine association are 
the high spin-down luminosity of the pulsar, the fact that the 
PWN’s cornet^ tailpoints back toward the SNR’s interior 
(seen in Figs. [2 and |2), and the edge-brightening and flatter 
spectral index of SNR G5.4-1.2 a long its western rim near the 
pulsar (Fig. m lFV^ et al. 11199 ' 4 ^ . In the absence of any con¬ 
straints on the PWN’s proper motion or the pulsar’s age, the 
simplest interpretation would seem to be that the pulsar was 
born reasonably close to the SNR’s geometric center, has then 
moved outward, overtaken the SNR, and in the process has 
re-energized the shell’s emissio n with its relativis tic outflow 
of particles and magnetic fields llShull et alJl989l) . However, 
in ^we have derived tp> 69-88 kyr if the pulsar was born 
near the SNR’s center. Since = 15 kyr ^ tp, the picture 
proposed above is problematic. 

In the following discussion, we consider three possible so¬ 
lutions'' to this difficulty; (1) the pulsar and the SNR are 
physically associated with an age tp « Tc, but the pulsar was 
born substantially offset from the SNR’s geometric center; (2) 
the pulsar and SNR are physically associated, but tp ^ r^; or 
(3) the two objects have no physical association. We conclude 
by considering the nature of G5.27-0.90, which we argue pro¬ 
vides additional clues to distinguish between these possibili¬ 
ties. 

3.1. Association Involving a Young Pulsar 

In the case of a physical association, but with a supernova 
explosion substantially offset from the geometric center of the 
SNR, our upper limit on tp can potentially be made consistent 
with the characteristic age. If the blast center of the SNR 
were considerably closer to the current location of the pulsar 
than the separation of 0 = 16.'l-20(6 assumed earlier, the 
predicted pulsar motion would be lower, and possibly could 
brought below the upper limit in Table [2 

A reasonable mechanism through which this large offset 
could have occurred is if the supernova occurred in the stellar 

*' We acknowledge the existence of a variety of other explanations that 
have been proposed_JorPSR_B1752_^4 and_^NRG5;4__L2j_™ch_as have 
been discussed bv IKundi l I I 992 11. IlstominI 119941) .I ^arsden et ^ 1 200 i ll and 
IShi & XiiliOOl . However, here we focus on three simple possibilities which 
can potentially be compared using available data. 


wind bubble of a moving progenitor star JGvaramadzel2002h . 
Expansion of a SNR into such a bubble could produce a rem¬ 
nant whose bl ast site is substantia lly separated from the geo¬ 
metric center JGvaramadzel2004l . In such a case, the pulsar 
could have been born inside the SNR but quite close to its 
present location, and the predicted proper motion would then 
be very low. 

This explanation was viable for the upper limits obtained 
by GFOO and TBG02. However, for our new, more constrain¬ 
ing proper motion measurement, the pulsar would need to 
have been born on the rim of or even outside the SNR to sat¬ 
isfy the requirement that tp k. Tc- While it is possible that 
the progenitor star was in the process of escaping its own 
wind bubble when the supernova occurred, such a system 
would produce an SNR significantly elongated and distorted 
along the east-west ax i s corresponding to the pulsy ’s motion 
llRoivczka et alJlT993t iBrighenti & D Ercolel 19941) . not con¬ 
sistent with the SNR morphology seen here. 

Furthermore, since the pulsar was kicked randomly away 
from its birthplace, a wide separation between the neutron star 
birth site and the SNR’s geometric center implies that there 
should most likely be a significant misalignment between 
the pulsar’s direction of motion (as inferred from the PWN’s 
morphology) and the ve ctor joining the SNR’s center to the 
pulsa r’s current position JGvaramadzel2002t iMigliazzo et al.l 
I 2 OO 2 I) . Specifically, for an explosion occurring near the edge 
of the cavity, there is a Ri 90% probability that the misalign¬ 
ment between the projections of these two vectors will be 
larger than the ss 20° observed. The mild misalignment be¬ 
tween the projections of these two vectors therefore makes a 
scenario in which the explosion was substantially off-center 
relatively unlikely. 

3.2. Association Involving an Old Pulsar 

As suggested by GFOO and lGvaramadzelll2004l) . the case for 
a pulsar / SNR association can be made if the characteristic 
age is a significant underestimate of the system’s true age. 
The solution proposed by GFOO predicts an evolved system 
consisting of an old SNR and a slow-moving pulsar, for which 
tp = 93 - 170 kyr and /r = 7- 10 mas yr“', consistent with the 
upper limits found here. However, we then need to explain 
why fp > Tc. 

To consider this possibility, we must reconsider the assump¬ 
tions underlying Equation Q. The braking index, n, is de¬ 
fined by the equation i> = -Ko", where o = 1/P. For spin- 
down via a magnetic dipole i n vacu o, we expect n = 3. Fol- 
lowing iBlandford & Romanii (Il988h . we separately define the 
“deceleration parameter”, h = vvliP-. If K and n are both 
constant, then h = n. 

There are two situations in which the characteristic age, = 
P/2P, can underestimate the true age as required here. The 
first possibility is that P ^ Pq and K is constant, but n = n<3. 
For PSR B1757-24, Equation 0 can only yield tp > 69 kyr 
for Po < 13 ms, and even then only for « « 1. In contrast, 
measureme nts of initial spin in most other systems s uggest 
Po > 20 ms JMigliazzo et alJ2002HKramer et alJ2003h . while 
in the few cases where n has been determined, it is found that 
1.4 < n < 2.9 (Livingstone et al. 2005, 2006, and references 
therein). It thus seems unlikely that a viable set of parameters 
can describe the system via Equation 0 . 

The alternative is that K is not consta nt. In this case. Equa¬ 
tion (p) no longer holds, and ii^n (see [Blandford & Romanii 
Il988^ . Either a changing magnetic held or a chang ing mo- 
ment of inertia can cause K to vary with time (see e.g.. lCamU(3 
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Il996l) . In this case, the pulsar spin-down is uncoupled from 
the star’s age, and rather traces the time scale on which K 
evolves. 

While it seems unlikely that the star’s moment of iner¬ 
tia could evolve substantially on this time-scale, the growth 
of the surface magnetic field is a possibility, as may re¬ 
sult from thermoelectric instabilities in the crust or diffu¬ 
sion of magnetic flux from the s tar’s interior (iBlandford et^ 
1983t iRuderman et al.l Il998t iKonenkov & Gennertl 1200 if 


Lin & Zhandl2004t . If we assume n = 3 but n < 3, then Tc 


underestimates th e true age, a nd the surface magnetic field, B, 
grows with time llEvnell20M) . Indeed ii < 3 is observed for 
all six pulsars in which this quantity has been accurately de¬ 
termined, while dB/dt > 0 is sugge sted in long-term timing 
signa t ures of seve ral other pulsars (iLvne et alJll996t ISmithl 
Il999t lLvnell2()04ll . While n has not been directly measured 
for PSR B1757-24, we can write as an order of magni¬ 
tude estimate that the time scale for field growth has been 
B/B ~ Tc « 15 kyr. 

We further note that this does not require us to invoke a 
unique effect to explain this specific pulsar / SNR associa¬ 
tion. Even though '^C ^ for many oth er pulsars in SNRs, 
these results are still co nsistent with « < 3 ilKasni et al 12001 lit 
iMigliazzo et alJl2002^ . Gradual magnetic field growth may 
thus be widespread, but is only noticeable in a system such as 
the Duck for which its comparatively large age has allowed 
this effect to accumulate. 

It is important to acknowledge that the large implied age 
{tp > 69 kyr) is beyond what is expected for the observ- 
able lifetime of a SNR, which is t ypically ~ 20-60 kyr 
llBraun et alJll989HErail et al.lll994al) . Thus even if magnetic 
field growth can explain the discrepancy between the pulsar’s 
characteristic age and that inferred from proper motion, the 
large implied age is also a potential issue for SNR G5.4-1.2. 
However, the faint eastern rim of the SNR is consistent with 
what is seen for older, undisturbed shells, and may not have 
been easily identified on its own if not for the brighter western 
side. Thus the large age required for the SNR is not inconsis¬ 
tent with the appearance of this half of the SNR. 

Eurthermore, if the pulsar and SNR are physically associ¬ 
ated, the discrepancy between the system’s inferred age and 
the expected lifetime for observable SNRs can be resolved 
by the argument tha t the pulsar is re-e nergizing the SNR as 
it passes through it llShull et al.ll1989ti . The re-energization 
hypothesis predicts that brighter emission with a flatter spec¬ 
trum should be seen along the western edge of SNR, and that 
these effects should peak where the pulsar would have crossed 
the rim. Indeed these phenomena are both observed (Eig. [0 
iBecker & Helfandlll985t lErail et aI1ll994bl) . However, a re¬ 
quirement of this theory is that synchrotron-emitting particles 
must diffuse sufficiently rapidly away from the pulsar around 
the rim of t he shell to produce the observ ed region of apparent 
interaction Jvan der Swaluw et alJ2nn2h . We can quantify this 
as follows. If T is the time elapsed since the pulsar first began 
to interact with the SNR shell and X is the distance which is 
traveled by particles along the shell’s circumference after in¬ 


jection by the pulsar wind, then we require X = (2KDuckT)*/^, 
where kducIc is the diffusion coefficient of synchrotron emit¬ 
ting particles as they move from the pulsar around the rim. 
Erom Eigure[2we estimate A « 10 pc and T ~ O.lf^ > 7 kyr. 
We thus require Kouck ~ 2 x 10^^ cm^ s“' to explain the SNR’s 
appearance. 

We can determine if this inferred diffusion coefficient is 


reasonable by considering the magnetic field configuration 
through which particles must propagate. If the SNR is 
in the radiative phase as argued by GEOO, then we expect 
strong compression at the SNR shock, and hence a mag¬ 
netic field which runs parallel to the shell. Indeed radio 
polarization measurements clearly demonstrate this geome¬ 
try to be present along the SNR’s western rim JMilne et alJ 
[1992). The mean free path in directions parallel to the mag¬ 
netic field lines is larger than the electron gyroradius by a 
factor pDuck = {SBo/Bo)~^, where Bq is the ambient magnetic 
field strength in the SNR sh ell and 6Bn is the amplitude of 
turbu lent fluctuations in Bq d.lokiDiill1987t lAchterberg et alJ 
l1994l) . We can thus write Kouck = iTDuckKB, where kb is 
the standard Bohm diffusion coefficient, with value kb « 
1 X cm^ s~* for the 327 -MHz image shown 

in Eigiire fTl ilvan der Swaluw et al.ll200^ . We then require 
X 10"^. In comnarison. lAchterherg et all 
(119941) find Pyoung ^ 2000 for young SNRs which have en¬ 
hanced turbulence, but pisM ^10^ for the ISM. If we assume 
a field strength Bq ~ 5 - 10 /rG, we thus find ryouck ~ Vism, 
consistent with SNR G5.4-1.2 being an old remnant with 
reduced turbulent amplitude which is merging into the ISM 
(see also lMoffett & Revnoldsl1994l) . We conclude that the re¬ 
energization hypothesis is consistent with the expected rapid 
diffusion of particles from the pulsar around the shell rim. 

As a final note, we point out that even if the pulsar is old, 
a small offset between the supernova explosion site and the 
SNR’s geometric center is still required. Specifically, the tra¬ 
jectory of the pulsar as inferred from the cometary tail of 
the PWN passes north of the SNR’s geometric center, indi¬ 
cating that the blast and g eometric centers d o not coincide 
(FigiS lErail et ai]ll994bh . lErail et akl lll994bl) proposed that 
a gradient in the density of the ambient interstellar medium 
(ISM) into which the SNR is expanding could produce an 
asymmetric expansion, indeed resulting in a small offset be- 
tween the explosion site and the SNR’s center. Alternatively, 
iGvaramadzel il20n4h has proposed that the misalignment be¬ 
tween the pulsar’s expected and inferred trajectories could be 
due to a progenitor star which explodes slightly offset from 
t he c enter of its wind-blown bubble (as already discussed in 
EH but for the case of a young pulsar with a subtantial offset 
from the bubble’s center). This again leads to the conclusion 
that the pulsar’s t rue age is considera bly in excess of its char¬ 
acteristic age ('see lGvaramad7el2004h . 

3.3. Chance Alignment 

The alternative simple explanation proposed by TBG02 is 
that SNR G5.4-1.2 and PSR B1757-24 have no physical con¬ 
nection. Given the relatively high density on the sky of both 
pulsars and SNRs in the inner Galactic plane, it is reasonable 
that tw o such objects could appear near each other in pro¬ 
jection dGaensler & Johnstonlll995h . Hi absorption is only 
able t o provide lower lim its on the distances to the SNR and 
PWN dPrail et al.l1994hl) . while the pulsar’s distance estimate 
as derived from its dispersion measure is model dependent 
and comes with significant associ ated uncertainties (see dis¬ 
cussion bv ICordes & Ea7i3l2002ft . Thus, although the dis¬ 
tances to the two objects are consistent, TBG02 point out that 
this is not necessarily strong evidence in favor of an associa¬ 
tion. 

Previous discussions had focused on the flatter spectrum 
and brighter emission of the SNR as evidence for an asso¬ 
ciation, as described above. However, TBG02 point out that 
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variations in spectral index and asymmetries in brightness are 
both common in SNRs. In particular, many SNRs are brighter 
on the side closest to the Galactic Plane, as observed here. 
While these are valid arguments, we note that in most other 
SNRs in which spectral index variations are observed, these 
changes are spread randomly o ver the SNR, rather than con - 
centrated in a particular region (lAnderson & Rudnicklll993l) . 
More specifically, the systematic trend toward flatter spectra 
as one gets closer to the pulsar position around the rim of 
G5.4-1.2 is difficult to explain if the SNR and pulsar are unas¬ 
sociated. 

3.4. The Nature of G5.27-0.90 

An important additional aspect of this discussion is the na¬ 
ture of G5.27-0.90. As TBG02 note, the very tight angular 
coincidence of PSR B1757-24 and G5.27-0.90, along with 
the cometary morphology of the latter, make it virtually cer¬ 
tain that at least these two objects are associated. 

It is important to bear in mind that Figure |2] shows only 
the westernmost extent of this source. The inset to Figure [I] 
demonstrates that eastward of the pulsar, this structure broad¬ 
ens into a larger nebula, approximately 100" x 100" in ex¬ 
tent. TBG02 interpret this overall morphology as a “Crab¬ 
like” PWN. In their model, the pulsar was born at the center 
of this larger nebula, and is now moving away from this site. 
However, there are difficulties with this interpretation. For the 
pulsar to have escaped from its own wind-driven bubble, the 
expansion speed of the PWN must have fallen well below the 
pulsar’s space velocity. This process can only occur in two 
scenarios. 

The first possibility is that as a PWN expands into the freely 
expanding ejecta of its associate d SNR, it will eventually col¬ 
lide with the SNR reverse shock (Rey nolds & Chevalieil19841 
Ivan der Swaluw e~ aillTiioIllBiondin et alJl2001l) . The result- 

ing interaction can reduce and even reverse the expansion of 
the PWN. Combined with the pulsar’s ballistic motion, this 
compression of the PWN produces a morphology in which 
the pulsar is at the tip of a trail of radio/X-ray emission, con- 
nected to a larger “relic PWN ”. This scenario was invoked by 
Ivan der Swaluw et alJ j2004l) to explain the morphologies of 
the PWNe seen in SNRs N157B and G327.7-1.1. The dif¬ 
ficulty with this interpretation here is that to have produced 
a reverse shock, the SNR needs to have interacted with and 
swept up a significant amount of interstellar gas. The fact that 
we do not see any radio emission from an associated SNR 
makes this interpretation problematic. 

The alternative possibility is that no significant outer blast 
wave was prod uced, as may be the case for the Crab N ebula 
and for 3C 58 llSlane et al.ll2004t ISeward et alJl200(Th . The 
PWN then interacts directly with the ISM, and has decelerated 
sufficiently as it sweeps up ambient gas that the pulsar is now 
able t o overtake it. In this case we can write dCastor et al.l 

[ml): 

( 2 ) 

where RpwN is the radius of the “relic” component of the neb¬ 
ula at time f, p is the ambient mass density, and where we 
have assumed that the pulsar blows a steady wind of luminos¬ 
ity E into a uniform surrounding medium. For t = Tc^ 15 kyr, 
E = 2.6 X 10^® ergs s“' and Rpwn ~ 1 pc, we then find an am¬ 
bient number density no ~ 500 cm“^. 

This high value can be ruled out by the morphology of the 
bow-shock component of the PWN. As discussed by GFOO, 
the stand-off distance between the pulsar position and the 


apex of the bow shock is r„ 0.04 pc. This position is set by 
pressure balance between the pulsar wind and the ram pres¬ 
sure of the surroundings, E /47rr^c « pV^, where V is the 
pulsar space velocity. The left hand term yields a pressure 
^ 5 X 10“'° ergs cm“°. For the density inferred above, we 
then find V ~ 10 km s“'. Not only is this velocity more than 
an order of magnitude slower than seen for young pulsars, but 
the solution is not self-consistent, since in 15 kyr, such a pul¬ 
sar would only have traveled ~ 0.1 -0.2 pc and so could not 
have moved outside its PWN bubble. 

We thus conclude that there is no obvious explanation for 
G5.27-0.90 if PSR B1757-24 is ~ 15 kyr old and is unasso¬ 
ciated with SNR G5.4-1.2. 

On the other hand, if there is a genuine association be¬ 
tween the pulsar and the SNR, then G5.27-0.90 could repre¬ 
sent the remnants of the interaction between the pulsar wind 
and the SNR shell. Ivan der Swaluw et al.l (l2003l) consider the 
interaction between a high-velocity pulsar and a SNR in the 
Sedov-Taylor phase of evolution. They show that as the pul¬ 
sar crosses the SNR, the drop in ram pressure from that in the 
SNR shell to that of the ambient ISM results in an expansion 
of the bow-shock structure. However, if G5.4-1.2 is in the 
pressure-driven “snowplow” phase of evolution (as argued_^ 
GFOO and as would be typical for the age inferred in TableQ, 
then the compression ratio of gas swept up by the forward 
shock will be much larger th an the standard facto r of four seen 
in the Sedov solution (e.g.. lBlondin et aPllOOSl) . Therefore, 
during the period in which the pulsar breaks through the rim 
of the shell, the ram pressure experienced by the bow shock 
drops dramatically, which in turn should produce a sudden, 
explosive expansion of the pulsar wind nebula. As the pulsar 
moves outward, the pulsar wind blows a “Crab-like” struc¬ 
ture into the ISM, which might possibly correspond to G5.27- 
0.90. Eventually the pulsar readjusts to its new, lower density 
environment, again forming a standard Mach cone, but con¬ 
nected to the larger relic structure left behind. This idea needs 
testing by full hydrodynamic modeling, but it provides a fea¬ 
sible explanation for the overall morphology of G5.27-0.90 in 
the scenario in which the pulsar is old and is associated with 
the SNR. _ _ 

We note that iGvaramadzel il2004l) has proposed an alternate 
explanation for G5.27-0.90, in which this source is a slowly 
expanding lobe produced by the collimated flow of hot gas 
which follows the pulsar as it punctures the SNR shell. This 
similarly requires a physical association between the pulsar 
and the SNR. 

4. CONCLUSIONS 

The upper limit on proper motion we find for the 
PWN G5.27-0.90 is the most restrictive one obtained to date 
for this system. We are able to reject the original expectation 
that PSR B1757-24 was born ~ 15 kyr ago and is associated 
with SNR G5.4-1.2, regardless of whether the corresponding 
supernova occurred at the geometric center of the SNR, or at 
a site substantially offset from this. 

Two possibilities remain. First, the pulsar and SNR are as¬ 
sociated and share an age > 70 kyr. In this case, the pulsar has 
a typical projected space velocity of < 330 km s“', which al¬ 
lows it to overtake its SNR at this stage and drive a bow shock 
through ambient gas. The sudden drop in pressure as the pul¬ 
sar crosses the SNR’s radiative shell might also explain the 
bulbous component of the PWN G5.27-0.90 seen between the 
pulsar and the SNR. The pulsar’s passage has re-energized the 
radio emission from the SNR through rapid diffusion of pul- 
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sar wind particles along magnetic field lines; without this new 
injection of particles, the SNR would be much dimmer, would 
have a steeper spectrum and generally would be more difficult 
to detect. The slight offset between the pulsar’s inferred tra¬ 
jectory and the SNR’s geometric center possibly results from 
a density gradient into which the SNR is expanding, or from 
a slightly offset explosion within a pre-existing cavity. The 
only difficulty with this picture is that the system’s age must 
then be many times larger than the pulsar’s characteristic age, 
which is not seen in other pulsar / SNR associations. The pro¬ 
posed explanation is that the surface magnetic field of the pul¬ 
sar is at the present epoch growing on a time scale of ~ 15 kyr. 
This effect is consistent with the properties of other pulsar / 
SNR associations, and is possibly also bein g seen in the long¬ 
term timin g signatures of several pulsars (iLvne et ^Il996t 
lEvnel2n04li . 

The alternative is that the pulsar is 15 kyr old as indicated 
by its characteristic age, and is unassociated with SNR G5.4- 
1.2. This explanation eliminates the need to invoke off-center 
cavity explosions, re-energized shells or growing magnetic 
fields to explain the observations. However, this interpretation 
offers no easy explanation for the morphology of G5.27-0.90 
which, contrary to the proposal of TBG02, cannot be easily 
explained as a relic nebula left behind at the pulsar’s birth site. 
This model also requires that the brightened emission and flat 
spectrum of the SNR near the pulsar be a coincidence. 

Erustratingly, a full understanding of the Duck remains 
elusive. However, future VLA or VLBA observations of 
PSR B1757-24 should eventually detect motion. Eurther- 


Achterberg, A., Blandford, R. D., & Reynolds, S. R 1994, A&A, 281, 220 
Anderson, M. C. & Rudnick, L. 1993, ApJ, 408, 514 
Arzoumanian, Z., Chernoff, D. F., & Cordes, J. M. 2002, ApJ, 568, 289 
Becker, R. H. & Helfand, D. J. 1985, Nature, 313, 115 
Blandford, R. D., Applegate, J. H., & Hemquist, L. 1983, MNRAS, 204, 1025 
Blandford, R. D. & Romani, R. W. 1988, MNRAS, 234, 57P 
Blondin, J. M., Chevalier, R. A., & Frierson, D. M. 2001, ApJ, 563, 806 
Blondin, J. M., Wright, E. B., Borkowski, K. J., & Reynolds, S. P. 1998, ApJ, 
500, 342 

Braun, R., Goss, W. M., & Lyne, A. G. 1989, ApJ, 340, 355 
Brighenti, F. & D’Ercole, A. 1994, MNRAS, 270, 65 
Camilo, F. 1996, in Pulsars: Problems and Progress, lAU Colloquium 160, 
ed. S. Johnston, M. A. Walker, & M. Bailes (San Francisco: Astronomical 
Society of the Pacific), 39 

Camilo, F, Bell, J. F, Manchester, R. N., Lyne, A. G., Possenti, A., Kramer, 
M., Kaspi, V. M., Stairs, I. H., D'Amico, N., Hobbs, G., Gotthelf, E. V., & 
Gaensler, B. M. 2001, ApJ, 557, L51 
Castor, J., McCray, R., & Weaver, R. 1975, ApJ, 200, L107 
Caswell, J. L., Kesteven, M. J., Komesaroff, M. M., Haynes, R. F, Milne, 
D. K., Stewart, R. T., & Wilson, S. G. 1987, MNRAS, 225, 329 
Cordes, J. M. & Lazio, T. J. W. 2002, astro-ph/0207156 
Fesen, R. A., Shull, J. M., & Saken, J. M. 1988, Nature, 334, 229 
Frail, D. A., Goss, W. M., & Whiteoak, J. B. Z. 1994a, ApJ, 437, 781 
Frail, D. A., Kassim, N. E., & Weiler, K. W. 1994b, AJ, 107, 1120 
Frail, D. A. & Kulkarni, S. R. 1991, Nature, 352, 785 
Gaensler, B. M. & Frail, D. A. 2000, Nature, 406, 158 (GFOO) 

Gaensler, B. M. & Johnston, S. 1995, MNRAS, 277, 1243 
Gvaramadze, V. V. 2002, in Neutron Stars in Supernova Remnants, ed. P. O. 
Slane & B. M. Gaensler (San Francisco: Astronomical Society of the 
Pacific), 23-26 

Gvaramadze, V. V. 2004, A&A, 415, 1073 
Hester, J. J. & Kulkarni, S. R. 1988, ApJ, 331, L121 

Hobbs, G., Lorimer, D. R., Lyne, A. G., & Kramer, M. 2005, MNRAS, 360, 
974 

Istomin, Y. N. 1994, A&A, 283, 85 
Jokipii, J. R. 1987, ApJ, 313, 842 

Kaspi, V. M., Gotthelf, E. V., Gaensler, B. M., & Lyutikov, M. 2001a, ApJ, 
562, LI63 

Kaspi, V. M., Roberts, M. S. E., Vasisht, G., Gotthelf, E. V., Pivovaroff, M., 
& Kawai, N. 2001b, ApJ, 560, 371 


more, an X-ray detection of G5.4-1.2 and of the eastern 
parts of G5.27-0.90 should provide additional co nstraints on 
the pro perties of the system (see discussion by iKasni et ^ 
|200^, while continued timing of the pulsar may be able to 
indicate the nature of the br aking torque and magnetic field 
evolution in this source (e.g. . ILvne et alJl996l) . 

Einally, other similar systems might provide additional 
clues. PSR B1951-I-32 is clearly i n the process of punctur¬ 
ing the shell of t he SNR CTB 80 (iHester & KulkarnilfTosl 
lEesen et alJl^ . while PSR J1016-5857 and SNR G284.3- 
1.8 have been prop osed as another such interacting system 
JCamilo et alj |2()0i[)^ O ver the ensemble of Galactic pulsars 
and SNRs. IShull et alJ ^1989^ predict approximately half a 
dozen systems in which the pulsar has recently penetrated the 
SNR shell. Searches for further interacting pairs may thus 
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TABLE 1 

Proper Motion, Velocity and Age Limits on G5.27-0.90 and 
B1757-24. 



GFOO 

TBG02 

This paper 

Proper motion (westward; mas yr“^) 

<24.8 

<37 

< 13.9 

Projected velocity (westward; km s“^)^ 

<590 

< 880 

< 340 

Age (kyr)’’ 

> 39-49 

>26-33 

>69-88 


Note. — All limits are given at 5-(T significance. 


^Assumes a distance of 5 kpc to the pulsar and PWN. 

^Assumes that the pulsar was bom © = 16f 1 -20f6 east of its current location, 
near the geometric center of SNR 05.4-1.2. 


Thorsett, S. E., Brisken, W. F., & Goss, W. M. 2002, ApJ, 573, LI 11 (TBG02) van der Swaluw, E., Achterberg, A., Gallant, Y. A., & Toth, G. 2001, A&A, 
van der Swaluw, E., Achterberg, A., & Gallant, Y. A. 2002, in Neutron Stars 380, 309 

in Supernova Remnants, ed. P. O. Slane & B. M. Gaensler (San Francisco: van der Swaluw, E., Downes, T. P, & Keegan, R. 2004, A&A, 420, 937 

Astronomical Society of the Pacific), 135-140 
van der Swaluw, E., Achterberg, A., Gallant, Y. A., Downes, T. P, & Keppens, 

R. 2003, A&A, 397,913 
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Fig. 1.— VLA observations of “the Duck”, with the SNR, PWN and pulsar indicated. The main panel shows a 327-MHz image of SNR G5.4-1.2 at a 
resolution of 55^^ x 41'', adapted from GFOO. Contours are at levels of 10, 25, 50 and 100 mJy beam“^, and the peak intensity is 150 mJy beam“^. The cross 
shows the approximate location of the center of SNR G5.4-1.2, while the horizontal line indicates the inferred direction of motion for the pulsar as implied by 
the morphology of G5.27-0.90. The inset shows a 8.5-GHz image of PWN G5.27-0.90, covering a 3!5 x 3!5 field at a resolution of 1" (image courtesy of 
NRAO/AUI and Dale A. Frail). 
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Fig. 2.— 8.5 GHz VLA observations of the western tip of PWN G5.27-0.90 at three epochs spread over 12 years. In each panel, the greyscale ranges from 
-0.1 to +0.3 mJy beam“^ (as shown in the wedge at right), while the contours are at levels of 0.6, 1.2, 1.8 and 2.4 mJy beam“^. Each image has been smoothed 
to a resolution of FWHM Of^85, as shown by the circle at the bottom right of each panel. The position of PSR B1757-24 at epoch 2002.2, as given by TBG02, is 
marked by an asterisk in the bottom panel, and has an uncertainty much smaller than the size of the symbol. The vertical dashed line shows the expected relative 
shift for westward proper motion of 24.8 mas yr“^ (corresponding to the 5-o- upper limit previously determined by GFOO). Note that the first epoch has poorer 
sensitivity than the second and third observations. 






